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ABSTRACT 





Satellite and surface data are utilized to analyze mesoscale 
and subsynoptic cloud and precipitation patterns. Digital GOES 
(Geostationary Operational Environmental Satellite) visible and 
infrared data are used to produce high resolution (4 n mi) 
satellite analyses of cloud amount, cloud type,  cloud-top 
temperature and height, and precipitation intensity for an 
approximate 1600 x 1600 n mi area over the northeastern United 
States and the western North Atlantic Ocean. Conventional 
surface observations, the ARS (Automated Radar Summary) chart and 
manual analysis of the imagery are used to evaluate the 
satellite-derived analyses for nine cases during the winter and 
spring 1985. 

The majority of cloud amount estimates for clear and overcast 
sky conditions are analyzed correctly; however, broken and 
scattered skies are underestimated. The general cloud type 
patterns are depicted accurately. More success occurs with 
uniform-textured clouds (i.e. nimbostratus, stratocumulus) and 
multi-layered clouds than with nonuniform-textured clouds (i.e. 
cirrus, cumulus). The majority of cloud-top temperature/height 
analyses are representative of the cloud types and patterns. 
Most precipitation areas are identified correctly; however, the 


intensity of the precipitation is underestimated. 
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Lo INTRODUCTION 


Satellite imagery and conventional surface data are primary 
data sources describing mesoscale weather features. Each provides 
a different perspective on mesoscale and subsynoptic processes. 
Continuous satellite imagery provides high resolution analyses of 
cloud systems from above while the coarser resolution 
conventional surface data describes the current ground-observed 
weather conditions. A basic problem is the integration of 
conventional surface data with satellite data. This problem is 
addressed here as both surface and satellite cloud data are 
utilized to analyze mesoscale and subsynoptic cloud and 
precipitation patterns. 

A cloud and precipitation analysis was developed (Wash et 
al., 1985) that uses digital GOES (Geostationary Operational 
Environmental Satellite) visible and infrared data to provide 
high resolution (4 n mi) satellite analyses. This analysis model 
uses Harris and Barrett's (1978) method to estimate cloud amount; 
a variation of Liljas' (1982) cloud-typing and precipitation 
intensity scheme; and Reynolds and Vonder Harr's (1977) 
bispectral method to derive cloud-top temperature. 

The first objective of this thesis is to attempt to 
generalize T Naval Postgraduate School (NPS) model (Wash et 
al., 1985) by modifying the infrared thresholds used in the 
algorithm to account for monthly variations in race and cloud- 
top temperature. The second objective is to monitor and evaluate 


the NPS-model's performance, using the newly-developed 


Deli 


thresholds, with data from the winter and spring seasons. A 
detailed evaluation of several summer cases was conducted by Wash 
et. al 1995) 

Chapter II presents recent cloud research involving digital 
satellite data. Chapter III describes the current NPS cloud and 
precipitation model and Chapter IV discusses how the algorithm 
used in the model was generalized for use in any season. The 
subjective and objective plan used in evaluating the  model's 
performance is described in Chapter V, with the results presented 
in Chapters VI, VII and VIII. Conclusions and recommendations for 


further investigation are found in chapter IX. 
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T RECENT CLOUD RESEARCH USING DIGITAL SATELLITE DATA 


During the past decade the determination of various cloud 
parameters using digital satellite data has been a vigorous area 
of research. Algorithms using bispectral thresholds, spatial 
coherence and radiative transfer computations have been developed 
to examine various cloud parameters such as cloud amount, cloud 
brightness and cloud-top temperature. 

Shenk, Holub and Neff 1(1975) developed a cloud type 
identification scheme for tropical ocean areas using Nimbus-3 
MRIR (Medium Resolution Infrared Radiometer) multispectral data 
from four spectral regions (.2-4 Um, 6.4-6.9 Um, 10-11 Um, 20-23 

um). Infrared thresholds used in making cloud type decisions were 
based ona radiative transfer model developed by Kunde (1967), 
which requires input of CEEE ne E CRE soundings and 
assumed maximum cloud-top heights for low (750 mb) and middle 
(450 mb) clouds in tropical ocean regions. It was assumed that 
cirrus and cumulonimbus cloud tops only existed above 300 and 250 
mb, respectively. The visible reflectance thresholds were based 
on the past experience of the authors; low visible thresholds for 
Clear ocean regions were based on an earlier scheme of Shenk and 
Salomonson (1972). The results of twenty case studies were 
successful; the derived cloud-type identification maps compared 
well to high resolution cloud photography. 

Minnis and Harrison (1984) examined the effects of extensive 

cloudiness over a large area using 8 km resolution GOES-E visible 


M55=-. 75m) and infrared (10.5-12.5 Um) data. They developed a 
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hybrid bispectral threshold method to extract the ftollowine 
radiative parameters: effective cloud-top temperatures and 
clear-sky temperature, cloud cover at three altitudes; cloud and 
clear-sky visible reflectance properties. Cloud amounts) and 
cloud-top temperatures are estimated for all three cloud levels 
(as well as total cloud cover) using a measured clear-sky 
temperature and infrared histogram containing infrared threshold 
counts corresponding to each level's cloud-base temperature 
determined by a lapse rate. The mean cloud brightness is 
computed using the infrared threshold-derived value of cloud 
amount, with the clear-sky brightness determined from minimum 
reflectance models and the measured average visible count. 
Surface-based cloud observations and photo analyses were used for 
evaluation purposes. Their results were reasonable;  satellite- 
derived cloud amount estimates were +.05 less than both ground 
observations and photo analyses. Agreement in estimating cloud 
cover was best at lower and higher cloud amounts and worst with 
cloud amounts near 50%. 

Coakley and Baldwin (1984) developed the spatial coherence 
approach for retrieving cloud cover for simple-layered systems 
using 8 km resolution NOAA-7 AVHRR (Advanced Very High Resolution 
Radiometer) GAC (Global Area Coverage) data. Their algorithm was 
designed to obtain cloud properties relating to climate studies 
rather than process imagery data. The spatial coherence method 
provides the fractional cloud amount and mean cloud-covered 
radiance. The plot of mean radiance versus standard deviation, 
which are calculated for each 2 x 2 array of GAC infrared data, 


has the shape of an arch. Uniform emission from all four pixels 
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is expected for either a cloud-free or completely cloud-covered 
oceanic region. Partially  cloud-covered pixel arrays rarely 
exhibit uniform emission. Therefore, the mean cloud-free and 
cloud-covered radiance values can be identified as the feet of 
the arch (low standard deviations) and the partially cloud- 
covered arrays make up the body of the arch. Coakley and Baldwin 
use only the 11 Um channel to calculate the means and standard 
deviations of the radiances to determine cloud cover but 
calculate the visible means and standard deviations to check the 
infrared values. The least reflective visible values should 
coincide with the warmest infrared values. The spatial coherence 
method is limited to ocean regions since land as a background is 
too variable. It is also restricted to simple-layered cloud 


systems and assumes that clouds are optically thick at the 


wavelength of observation. Therefore, the spatial coherence 
method is inapplicable to baroclinic zones (too complex), ITCZ 
(non-layered structure) and cirrus clouds (too transparent). For 


a three-day period over a 250 km area, Coakley and Baldwin found 

the typical uncertainty in estimating cloud cover to be 14%. 
Platt (1983) examined the bispectral method for determining 

cloud parameters as well as separating broken and semitransparent 


cloud types using GMS-1 visible (.5-.7 Um) and infrared (10-12pm) 


data. Bispectral data are used to estimate cloud cover and 
optical depth, and to identify different cloud types. Unbroken 
clouds with variable optical depths (i.e. cirrus) are 


distinguished from broken clouds with uniform optical depths 


(i.e. stratocumulus) by comparing the shapes of bispectral curves 
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(infrared brightness temperature vs. visible albedo). Platt 
discovered the following relationships: albedo increases rapidly 
with decreasing infrared brightness temperature when the optical 
depth approaches unity; the albedo and infrared brightness 
temperature increase with cloud amount. Platt concluded that the 
best three-dimensional cloud representation would be a 
combination of the bispectral and spatial coherence methods. 


Tsonis (1984) developed an analysis scheme which identifies 


general sky and surface conditions. His classification groups 
are clear skies/snow cover, clear skies/no snow cover and 
high/low broken clouds and overcast. He uses GOES visible (.7 


Hm) and infrared (11.5 um) data with spatial resolutions of 4 x 
4 km and 8x 8 km, respectively. The infrared raw counts are 
converted to temperatures and the visible raw counts are 
normalized to account for sun angle. The two primary procedures 
in his separation method are segmentation and classification. 
Tsonis segments the image into spatially continuous clusters and 
classifies these groups based upon spectral thresholds and 
spatial textural features. He defined the cloud/no cloud 
threshold range to be between 24 and 27 raw GOES sensor counts 


(corresponding to albedos of 0.12 ~ 0.19) depending upon the 


surface characteristics. Tsonis verified his classification 
results with ground synoptic reports. His agreement percentage 
with ground truth data was 72%. This percentage increased to 


approximately 87% when all clouds were considered as one 


category. 
The algorithm used in the NPS model is most closely related 


to the method of Tsonis (1984). Both use spectral thresholds and 
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spatial textual characteristics of the visual data to determine 
various cloud parameters and the satellite-derived analyses are 
verified with ground synoptic observations. Like Shenk et. al. 
(1975) and Platt (1983) a cloud type identification scheme based 
on infrared and visible thresholds is used in the NPS model. 
Coakley and Baldwin's (1984) spatial coherence scheme and Minnis 
and Harrison's (1984) hybrid bispectral threshold method were 
designed for climatological studies; the NPS model, which is 
described in the following chapter, was developed to analyze 


mesoscale real-time imagery. 


17 


TIS THE GEOSTATIONARY SATELLITE CLOUD AND PRECTP TRA aa 
ANALYSIS MODEL 





A general description of the automated NPS cloud and 
precipitation model is presented. (For further information, refer 
to the technical report of Wash et al. (1985) or Moren (1984).) 

The NPS model uses operational digital satellite data from 
the visual and infrared channels from the GOES Visual-Infrared 
Spin Scanner Radiometer (VISSR). The NPS model creates real-time 
analyses of cloud amount, cloud type, cloud-top temperature, 
cloud-top height and precipitation intensity. The model 
processes the satellite data, performs statistical calculations, 
and produces cloud and precipitation analyses. 

Visual raw data counts are converted to albedos using a 
normalization scheme that corrects for sun angle and anisotropy 
(Muench and Keegan, 1979). Infrared data counts are converted to 


temperatures using the GOES sensor conversion table (Corbell et 


al., 1978). For each 2 x 2 grid of visual data, an average 
brightness and standard deviation is computed. The average 


Visible brightness and corresponding infrared temperature are 
used in the final cloud and precipitation analyses. 

The hybrid bispectral threshold method of Reynolds and 
Vonder Haar (1977) is used to estimate cloud amount. cloud 
amount is determined by comparing the visual albedos with a pre- 
established cloud/no cloud threshold of 0.15. If the albedo is 
greater than or equal to the 0.15 threshold, the pixel is 
considered cloudy. The cloud/no cloud threshold was lowered from 


0.17, the value used in the summer evaluation (Wash et__al,, 
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1985), to 0.15 so that more small cumulus could be identified. 
The 0.15 threshold compares well with the 0.13 threshold used by 
Keegan and Niedzielski (1981), who looked at 1977-78 
autumn/spring/summer data over the northeastern United States, 
and the 0.12 - 0.19 albedo range used by Tsonis (1984). An 
average cloud amount for each 2 x 2 grid is calculated (Harris 
and Barrett, 1978; Fye, 1978). A larger analysis array would 
provide a more complete range of cloud amounts but take longer to 
execute. 

Cloud type is determined through a series of threshold tests 


that use both infrared temperatures and visual albedos (Liljas, 


1982). The infrared thresholds have been modified to compensate 
for the monthly variation in the infrared temperatures. These 
modifications are discussed in the following chapter. In 


addition, a texture test that computes the standard deviation of 
each 2 x 2 grid of visual data is performed to differentiate 
between cumuliform and stratiform clouds which often have similar 
visible and infrared (11 um) radiative properties (Liljas, 1984). 

The precipitation estimation method, a modified version of 
Liljas' (1982) threshold technique adopted from the results of 
Muench and Keegan (1979), applies the cloud type information to 
the precipitation areas. The precipitation module is activated 
if the cloud-type decision is nimbostratus or cumulonimbus. The 
magnitude of intensity is dependent upon the infrared and visual 
values; the colder infrared temperatures and brighter visual 
albedos produce heavier intensities (Liljas, 1984). The infrared 


precipitation/no precipitation threshold, which varies from 247.5 
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K to 254 K, compares well with the 246 K to 253 K temperature 
range of Paul (1983). 

Cloud-top temperatures are estimated using a bispectral 
approach (Reynolds and Vonder Haar, 1977) that provides a better 
analysis of cloud-top temperatures for partially-filled fields of 
view and cirrus layers since the amount of surface radiance is 
estimated. Cloud-top temperatures are derived from the computed 
radiance, which takes into account surface and cloud radiance, 
emissivity and cloud amount. Finally, the cloud-top temperatures 
are equated to the pressure level heights from the Fleet 
Numerical Oceanography Center (FNOC) grid-point temperature 
profiles. 

The following chapter describes how the model is made 
applicable for use in various seasons. The infrared thresholds 
used in the cloud-typing scheme were modified to compensate for 


the monthly variation in the infrared temperatures. 
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DV. GENERALIZATION OF ALGORITHM 


The more difficult objective of this research effort is to 
generalize the algorithm used in the NPS model so that it can be 
applied for any month. It was observed from the analysis of the 
winter and spring cases, discussed in the following chapters, 
that the cloud infrared temperatures become progressively warmer 
from January to May. The infrared thresholds specifying various 
cloud groups need to be modified to compensate for this warming. 
Manual nephanalysis and conventional surface data were used to 
tune the initial winter and spring case studies to arrive at a 
set of infrared thresholds for each case. Independent cases are 
used to test the new infrared thresholds. 

The visible thresholds as well as the infrared thresholds 
used in this investigation (Figs. 1-5) differ from the earlier 
summer values (Fig. 6). The following modifications were made to 
the visible thresholds: The cloud/no cloud threshold was lowered 
from 0,17 to 0.15; the visible threshold differentiating cirrus 
from nimbostratus/multi-layered or altostratus/cumulus congestus 
was lowered from 0.55 to 0.32; the upper visible threshold 
defining stratus was raised from 0.55 to 1.0 while the lower 
visible threshold defining stratocumulus was lowered from 0.55 to 
0.15 (stratocumulus and stratus are now only distinguished in the 
infrared spectrum); nimbostratus/multi-layered and cumulonimbus 
are now differentiated by a visible threshold 0.88 (> 0.88 
results in a cumulonimbus classification) but can have similar 


infrared characteristics. These modifications resulted from the 
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investigation of the analysis errors that occurred in the summer 
cases discussed in Wash et al. (1985). The visible thresholds 
are not altered throughout this study. 

The infrared thresholds defining low-level cloud types 
stratocumulus, cumulus and stratus exhibit the most variation 
from Monthy to mondi Liljas (1984) found that stratus cloud-top 
temperatures may vary up to twelve degrees during the summer. 
Maturi and Holmes (1985), who developed monthly infrared imagery 
enhancement curves for sea fog, used enhancement curves that 
varied from 276 K - 281 K in April to 287.5 K - 293 K in August. 
A threshold scheme, which utilizes real-time satellite data, 
takes into account daily temperature fluctuations; an infrared 
enhancement method, which uses pre-determined enhancement curves, 
may not identify small temperature pertubations. 

In the 17 Jan study, the infrared temperatures which classify 
stratocumulus are 256-266 K (Fig. 1). In the 28 Mar case, 
stratocumulus is classified between 286-265 K (Fig. 3) which is 
at least nine degrees warmer than the stratocumulus in the 17 Jan 
case. This warming trend is evident in the scattergram (Fig. 7) 
in which the visible and infrared raw data counts (lower infrared 
count for warmer temperature) corresponding to surface reports of 
stratocumulus are plotted for each case. The stratocumulus of 
the 28 Mar case (represented by the symbol X) 1S warmer than the 
stratocumulus of the 17 Jan case (represented by the symbol 0). 
A similar temperature variation occurs with stratus and cumulus 
clouds. The stratus and cumulus clouds appear to be at least ten 
degrees warmer in May than January; this is illustrated in Fig. 


DI 
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A dramatic increase in infrared r TE from January to 
May does not occur with middle-to-upper level clouds. The raw 
infrared counts of the multi-layered reports for each case 
exhibit a random pattern (Fig. 9) rather than a January to May 
warming trend. Therefore, the 247.5 K temperature threshold 
defining nimbostratus/multi-layered and cumulonimbus in the 17 
JAN case is only raised by 6.5 degrees to 254 K for the 
subsequent cases. The small threshold adjustment also applies 
for cirrus and altostratus/cumulus congestus. 

The significant temperature difference between months found 
with low-level clouds is due to the increase in surface 
temperatures from January to May. This surface and low cloud 
warming trend is illustrated in the three-dimensional cluster 
diagrams (Figs. 10-12). The number in the darkened regions 
represents the average frequency of occurrence (higher frequency 
in the darker regions) of each infrared and visible (scaled to 
the infrared data by multiplying by four) count combination that 
occurs within the image. For example, inthe January cluster 
diagram (Fig. 10), there are approximately 157 occurrences of an 
infrared count of 75 (293 K) and visible count of 60 (15 x 4; 
albedo of 0.02). The blank horizontal lines represent raw counts 
that are missing in the infrared imagery. The cluster diagrams 
separate the satellite digital data into clusters of similar 
visible and infrared raw counts. Each cluster can be identified 
in the images. For example, the cluster in the middle of Fig. 10 
(17 Jan) can be identified as stratocumulus clouds since the 


stratocumulus observed in the GOES imagery (Figs. 13 and 14) 
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possess these visible and infrared raw count values. The cluster 
diagrams indicate that upper-level clouds do not have much 
monthly temperature variation. A similar clustering scheme used 
by Hawkins (1980) separates only infrared imagery into different 
radiative regions in which attributes are evaluated for the 
various clusters. 

The performance of the NPS model, using the new infrared 
thresholds, was monitored and evaluated for several 1985 winter 
and spring case studies. The evaluation plan and results are 


described in the following chapters. 
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V. EVALUATION PLAN 


The eastern United States served as the evaluation region for 
this investigation. This geographical location was selected due 
to abundant verification data, variety of meteorological 
phenomena (including coastal, land and cloud features) and the 
ability for direct comparison to the summer results for the same 
area (Wash et al., 1985). 

GOES visible and infrared data, which are collected at the 
Naval Environmental Prediction Research Facility (NEPRF) in 
Monterey, California, are modified to provide a center point at 
35° N 80°W for a 512 x 512 array with 2 x 2 n mi visual (2 x 4 n 
mi infrared) resolution. The 512 x 512 array is divided into 
sixteen 64 x 64 grids. For each of the sixteen center points of 
the grids, a 1200 GMT surface and upper-level temperature profile 
is obtained from the Fleet Numerical Oceanography Center (FNOC). 
The final output cloud and precipitation fields are 256 x 256 
grid arrays which approximately covers a 1024 x 1024 n mi area. 

The verification data network over the eastern United States 
consists of National Meteorological Center (NMC) facsimile 
charts, FNOC synoptic land, ship and surface airway reports, and 
manual nephanalysis. 

Nine winter and spring cases were chosen to monitor and 
evaluate the performance of the NPS cloud and precipitation model 
using the new infrared thresholds. Five cases were used to 
create the new winter and spring thresholds used in the algorithm 


and four were used as independent test cases. Data collection 
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consisted of capturing coincident infrared and visible data as 
well as the verification data for 1800 GMT 17 Jan, 20 Mar, 28 
Mar, ll Apr, 13 May (initial cases) and 20 Feb, 27 Mar, 16 Apr, 
15 May, 1985 (independent test cases). A subjective and 
statistical evaluation was performed on each case in which the 
NPS model output was compared to the conventional data. 

Direct comparison between ground-observed sky conditions and 
satellite digital data is a difficult task. Tsonis (1984) cites 
the following difficulties: 


An observer reports the sky conditions for an area that 
is obviously larger than the area represented by 
point...Therefore, comparison between the satellite 
classification (corresponding to a point that coincides 
with the geographical location of the synoptic station) 
and the station's report may be quite 
inadequate...Furthermore, any classification by an 
Observer, as "higher" or "lower" clouds is less objective 
than that by a satellite...The following reasons are 
given:the observer...looks at a much larger area, seeing the 
clouds in that area at different angles. Individual 
reactions, as well as human error, may contribute to the 
classification. The satellite, on the other hand, will 
indicate higher or lower clouds by its infrared count al 
each point;such a procedure is very objective. Another fact 
that may give rise to problems is that the observer looks at 
the clouds from below, whereas the satellite looks from 
above. 


These ground-satellite perspective difficulties encountered in 
this investigation will be discussed further in the following 


chapters. 
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VI. SUBJECTIVE EVALUATION RESULTS 


The subjective evaluation results of the model's performance 
in the 17 Jan case are presented in full along with an overview 
of its performance on the other days studied. In general, the 
model produced accurate cloud amount, cloud-type, cloud-top 
temperature/height and precipitation intensity maps in each case. 

In this January case, the primary synoptic features are a 
warm front, which extends from a 1000 mb low pressure center 
located at 379N 750W, and a cold front that extends from Florida 
to the Gulf of Mexico. Extensive cloudiness and precipitation 
are associated with this frontal system as evident in the GOES 
visible and infrared imagery (Figs. 13 and 14). A secondary 
feature is an occluded front, extending into eastern Pennsylvania 
associated with a 996 mb low pressure center, which produced snow 
and snow showers throughout the region. 

Each cloud amount category is depicted by a particular color. 
In this case, the NPS model depicts the extensive overcast (blue) 
covering the majority of the region as well as the broad clear 
zone (gray) south of the front over the southern tip of Florida 
and adjacent waters (Fig. 15). Broken/overcast (dark green), 
broken (light green) and scattered (red) sky conditions are 
underestimated in the southeast quadrant. 

The estimation of cloud amount is directly related to the 
Size of the array used in making the cloud amount decision, which 
is determined by the number of pixels with albedos greater than 


the 0.15 cloud/no cloud threshold. A 2 x 2 matrix (approximately 
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4 n mi) of visible data used in the model allows only four 
fractional cloud amount estimates: 25 50%, 75%, 100%. 
Increasing the number of albedos (larger field of view) used in 
making the cloud amount decision would yield a better range of 
fractional cloud amounts which would lead to more estimates of 
broken and scattered skies. 

The eight cloud type classifications of the 17 Jan case, 
which are represented by a different color, are presented in Fig. 
16% The extensive stratocumulus field (light green) associated 
with the cold air to the rear of the baroclinic zone is correctly 
depicted by the model. The cirrus (red) and altostratus (orange) 
located in the northern quadrant of the image are also classified 
correctiy. The model's analysis of the large stratus (yellow) 
and cumulus field (dark green) located in the southwest quadrant 
is verified by manual imagery analysis. The NPS model also 
indicates correctly the line of cumulonimbus (dark blue) over 
northern Florida. 

The model produced a different cloud-type analysis than the 
ground observations in the eastern quadrant due to the difference 
in perspective between the ground and the satellite. The model 
cloud-typing decision for this area is nimbostratus/multiple 
layers (blue) while many surface stations report stratocumulus. 
The infrared temperatures range from 251-235 K and the albedo 
values are all greater than 0.55. The satellite senses cold 
infrared temperatures, indicating middle/high clouds, while the 
surface observer views only the low cloud layer. The ground 


observer can not distinguish the upper cloud layers. 
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The NPS cloud-top temperature analysis (Fig. 17) depicts a 
range of temperatures (220-290 K) which corresponds to the wide 
variety of cloud layers and types found in this 17 Jan case. 
Each of the nine colors corresponds to a 10 K cloud-top 
temperature interval which ranges from 210-219 K (dark blue); to 
290-300 K (gray). Over northern Florida the analyzed cloud-top 
temperatures are 220 K (blue) and 230 K (light blue), which 
verifies with the surface reports of cumulonimbus and the 1835 
GMT Automated Radar Summary (ARS) chart (Fig. 18). 

The cloud-top height analysis (Fig. 19) follows the 
temperature analysis since cloud-top height is a function of 
cloud-top temperature. Each of the nine colors corresponds to a 
100 mb cloud-top height interval which ranges from 100-199 mb 
(dark blue) to 900-999 mb (dark red). 

The NPS model defines two precipitation areas: northern 
Florida and along the eastern seaboard. Three colors are used to 
describe precipitation intensity: red, green and blue denote 
light, moderate and heavy intensity, respectively. The model 
distinguishes the squall line over Florida with estimates of all 
three degrees of intensities (Fig. 20 ). This is verified by the 
1835 GMT ARS chart (Fig. 18) where light rain (red) to heavy 
thunderstorms (blue) are detected by the radar. Light (red) and 
moderate (green) precipitation is associated with the frontal 
zone. This is also verified by the 1835 GMT ARS chart but only 
two coastal surface stations (indicated by the arrows) report 
precipitation in this area (Fig. 21). The model overestimates 
the frontal precipitation. Overall, the NPS cloud and 


precipitation model performed well in the 17 Jan case. 
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The spring cases (20 Mar, 28 Mar, 11 Apr and 13 May) were 
also evaluated subjectively. The highlights of the evaluation 
are presented for each cloud parameter: cloud amount, cloud type, 
cloud-top temperature/height and precipitation intensity. 

Cloud amount estimates are generally accurate for overcast 
and clear skies. The broken/overcast, broken and scattered 
situations tend to be underestimated since only four visible 
albedos are used in determining cloud amount. 

Some small cumulus elements are not identified as clouds 
since their albedos are below the 0.15 visual cloud/no cloud 
threshold. The small cloud elements do not fill the sensor's 
field of view; therefore, surface reflectivity influences the 
cloud/no cloud decision and cumulus is classified as clear skies. 
This is evident in the 20 Mar case (Fig. 22) where some small 
cumulus in Florida, as seen in the GOES visible and infrared 
imagery (Figs. 23 and 24), are not identified. 

Cloud amount is overestimated when snow Cover is mistaken for 
cloud cover (Minnis and Harrison, 1984). An example of this 
situation is found over Ontario north of Lake Erie in the 11 Apr 
case. This area is classified as cloudy (Fig. 25) due to bright 
albedos (> 0.15) and cold infrared temperatures (270 K) but is 
really clear and snow-covered as seen in the GOES visible and 
infrared imagery (Figs. 26 and 27); this is verified by the snow 
cover map UD TELAS. The AFGWC Automated Cloud Analysis Model 
does not calculate cloud amount whenever snow is present at a 
grid point; clouds can not be discerned from snow (Fye, 1978). 


Generally, the cloud/no cloud decision is made correctly in the 
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spring cases but the scattered and broken classifications are 
underestimated. 

In the spring studies the cloud type classifications that are 
the most successful are large areas of uniform-textured clouds 
such as nimbostratus and stratocumulus. For example, in the 28 
Mar case (Fig. 29) the model identifies the multiple-layered 
clouds (blue) of the baroclinic zone, and the stratocumulus 
(light green) and cumulus (dark green) to the rear of the front, 
as evident in the GOES visible and infrared imagery (Figs. 30 and 
51). 

Common cloud-typing errors are: (l) clear skies with snow- 
covered ground classified as stratocumulus and cumulus clouds; 
(2) semitransparent cirrus classified as clear skies or low 
cloud; (3) cumulus classified as clear skies. 

In the ll Apr case the NPS model specifies the clear snow- 
covered areas in the northwest quadrant as stratocumulus (light 
green) (Fig. 32) since both the visible threshold (0.15) and the 
infrared threshold (291 K) are exceeded. 

Another common misclassification occurs with cirrus-type 
clouds. Semitransparent cirrus, which is located in the 
Carolinas and Tennessee, is identified as clear skies in the 11 
Apr case (Fig. 32). The surface radiation reaching the sensor 
results in extraordinarily warm infrared temperatures which 
causes the misclassification. 

Small cumulus are sometimes classified as clear skies. This 
is illustrated in the 13 May and 28 Mar cases. In the 13 May 
case the small cumulus in Florida and the southwest quadrant 


(evident in the GOES visible and infrared imagery, Figs. 33 and 


aya 


34) are not identified by the model (Fig. 35). In the 28 Mar 
case the Cumulus along the Gulf Coast (evident in the GOES 
visible and infrared imagery, Figs. 30 and 31) are classified as 
clear skies by the model (Fig. 29). The small cloud elements do 
not fill the sensor's field of view which allows surface effects 
(warm infrared temperatures and dim albedos) to influence the 
cloud-type identification decision. 

The cloud-top temperature/height analyses are fairly accurate 
in the spring cases. For example, in the 28 Mar case the cloud- 
top temperature analysis (Fig. 36) depicts the various cloud 
types associated with the baroclinic zone. Cloud-top temperatures 
of 220 K (blue) and 210 K (dark blue) represent cirrus-type 
clouds. Temperatures ranging between 230 K (light blue) and 260 
K (light green) indicate multiple-layered clouds and the 
stratocumulus and cumulus to the rear of the system have cloud- 
top temperatures between 270 K (yellow) and 280 K (red). The 
only cloud-top temperature errors occur with semitransparent 
cirrus and small cumulus for similar reasons as discussed above 
(i.e. the influence of surface radiation). 

Precipitation delineation is reasonable for the majority of 
spring cases. For example, in the 20 Mar case the large 
precipitation area identified by the model in the central and 
southwest quadrant (Fig. 37) is verified by the 1835 GMT ARS 
chart (Fig. 38) and surface observations. The small precipitation 
region over Maryland is also correctly identified. 

A common precipitation error throughout the investigation is 


also illustrated in the OMA The radar detects 


32 


precipitation occurring over western Pennsylvania but the model 
misses this area since the infrared temperatures are below the 
precipitation/no precipitation threshold (Fig. 39). 

Precipitation intensity is more difficult for the model to 
estimate and to be evaluated by an investigator. Lovejoy and 
Austin (1979) conclude that GOES infrared and visible data are 
good for determining precipitation areas but poor for determining 
intensities. The biggest error in precipitation intensity is 
found in the 11 Apr case. The NPS analysis (Fig. 40) indicates 
light precipitation (red) and the radar detects moderate rain and 
rain showers off the northeast coast (Fig. 41). Overall, the NPS 
model can determine precipitation distribution and does a fair 
job in estimating intensity. 

The subjective evaluation results of the model's performance 
with winter and spring data were described in this chapter. The 
objective evaluation in the following section presents additional 


outcomes and supports the preceding conclusions. 
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VIT. OBJECTIVE EVALUATION RESULTS 


A statistical evaluation is performed on the cloud type and 
precipitation parameters for each winter and spring study. The 
NPS model cloud type and precipitation analyses are compared to 
conventional surface data. 

The verification network used in the evaluation is a 
combination of 1800 GMT synoptic land, airways and ship data for 
an area encompassing the satellite image, 25-45° N and 90-70? w. 
Tsonis (1984) also uses ground synoptic stations to verify his 
satellite classification scheme. 

The surface data processor used, an updated version of the 
scheme developed by Wyse (1984), only processes cloud type and 
precipitation information. The surface data analysis methods 
utilized in the processor are similar to those of Hahn et al. 

(1984) and are summarized as a flow chart in Fig. 42. The 
surface observations of sky cover, present weather and low, 
middle and high cloud type determine the cloud-typing decision. 
Hahn et al. (1984) and the NPS processor discard the surface 
report if the sky cover information is missing, if the sky is 
obscured and there is no precipitation reported, or if there is 
sky cover reported but the cloud type information is either 
missing or zero and no report of weather. The NPS version 
classifies those reports with sky cover less than or equal to 
three tenths as clear skies. Tsonis (1984) uses this same 
criterion in his classification scheme. If just the cloud 


information is missing but present weather exists, the cloud type 


34 


is determined by the weather parameter. For example, if the 
cloud type group is missing but rain is reported, the cloud type 
decision is nimbostratus. 

Synoptic and airways observations are classified into eight 
categories (Ci, As, St, Sc, Cu, CC, Ns and Cb) corresponding to 
the World Meteorological Organization (WMO) code (Tables 1-3). 
The decision process is summarized by Fig. 43. If two or more 
cloud layers are observed, the report is classified as multiple 
layers (ML). Two differences between the actual NPS cloud-typing 
decision process and the Hahn et al. (1984) scheme are that the 
NPS processor distinguishes between stratus and stratocumulus and 
cumulus and cumulus congestus. 

As seen in Tables 1-3, some cloud type categories are 
misleading; therefore, special processing decisions are needed. 
For example, low cloud type 7 ( stratus/cumulus of bad weather 
or both usually below altostratus or nimbostratus) incorporates a 
middle-level cloud (altostratus, nimbostratus). Therefore, if the 
cloud cover exceeds five tenths, the classification is multi- 
layered clouds; if it is less than or equal to five tenths, the 
resulting cloud classification is stratocumulus. Low cloud type 
8 (cumulus and stratocumulus) does not distinguish between the 


two categories; therefore, these low cloud type reports are 


ignored. Another example occurs with middle cloud type 7 
(altocumulus, etc.). Altocumulus is not distinguished from 
nimbostratus; therefore, the cloud decision is 


altostratus/nimbostratus when no precipitation/precipitation is 


reported. 
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The precipitation decision is based on the present 
weather parameter. If the present weather is a WMO code between 
50 and 99, precipitation is occurring at the time of the 
observation. 

Each valid surface observation of cloud type is compared to 
the corresponding NPS model cloud type decision as well as the 
predominant cloud type of a 7 X 7 matrix which covers an 
approximately 28 x 28 n mi area. A voting procedure is used on 
the matrix to determine the cloud type that occurred most 
frequently in the satellite matrix. A straight percentage is 
computed for the number of agreements between the surface 
observation and the satellite information (assuming the surface 
report to be correct): PERCENTAGE CORRECT = NUMBER OF AGREEMENTS 
/ TOTAL NUMBER. A percentage correct is computed for clear 
conditions, eight individual cloud type categories (Ci, As, St, 
Sc, Cu, CC, Ns and Cb) and three cloud group categories (low, 
middle and high). Surface reports of multi-layered clouds are 
evaluated ina subjective rather than an objective manner since 
the NPS cloud-typing routine does not separate multi-layered and 
nimbostratus clouds. 

The multi-layered cloud type category comparison is made 
between the 7 x 7 satellite matrix and the corresponding surface 
observation. Surface reports of multiple layers are placed into 
five separate categories according to the thickness properties of 
the particular cloud types. 

Reports of dense cloud types at two or more levels are 
classified as multi-layered (i.e. stratocumulus, WMO low cloud 


type 5 and dense cirrus, WMO high cloud type 2). Ground 
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observations of semitransparent cloud types at all three levels 
are categorized as thin clouds. Semitransparent middle and/or 
lower level cloud types, but dense clouds at upper levels are 
considered high clouds and thin upper and/or lower layer cloud 
types, but opaque middle cloud types are designated as middle 
clouds. Finally, semitransparent cloud types at upper and/or 
middle levels, but thick at the lower level are classified as low 
clouds. 

Each 7 x 7 matrix corresponding to a surface report of 
multiple cloud types is also categorized. If two or more cloud 
types at different levels occur in the matrix, the satellite 
cloud-typing decision is multi-layered. If the matrix contains 
mostly zero values (clear conditions), the resulting 
classification is thin. If the majority of cloud types are 
cirrus, the decision is high cloud. The matrix is specified as 
middle cloud if the majority of cloud types are middle cloud 
(altostratus,  nimbostratus) and it is classified as low cloud if 
the predominant cloud types are low clouds (stratocumulus, 
stratus, cumulus, cumulus congestus). 

Assuming the surface report of multi-layered clouds to be 
correct, each satellite matrix is compared to the ground truth 
data. If the surface and satellite do not have the same multi- 
layered properties, (i.e. surface report is classsified as thin; 
satellite is specified as multi-layered) the pair is labelled 
mismatch. 

For each case, precipitation occurrence is also evaluated 


Statistically. The surface report of precipitation/no 
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precipitation is compared to the result of the voting scheme 
performed on the 7 x 7 matrix in which twenty-four precipitation 
votes out of forty-nine pixels constitutes precipitation. 

Each surface report of precipitation is compared to the NPS 
model's precipitation decision and a percentage correct is 
computed: PERCENTAGE CORRECT FOR PRECIPITATION REPORTS = NUMBER 
OF AGREEMENTS BETWEEN SATELLITE AND SURFACE / TOTAL SURFACE 
REPORTS OF PRECIPITATION. The same evaluation procedure is used 
for surface stations where no precipitation is occurring. 

The results of the cloud type and precipitation evaluation 
for the winter and spring cases are found in Tables 4-7. Table 4 
contains the statistical results of the specific cloud type 
categories, the group cloud type statistics are found in Table 5 
and the multi-layered results are presented in Table om 
Precipitation results are presented in Table 7. 

The most successful cloud type classification is the clear 
skies category. The percentage correct (one to one 
correspondence between the NPS model output and surface report) 
ranges from 71% (11 Apr) to 1003 (17 Jan) and a total value of 
91% for all cases. In each clear sky case, the 7 x 7 vote 
produces as good or better results than the single pixel 
comparison (Tables 4 and 5). 

Reasonable results are obtained for the stratiform cloud 
types, stratocumulus and nimbostratus, for both direct and matrix 
comparison methods. In most instances these cloud types are 
specified correct nearly fifty percent of the time. These 


uniform-textured clouds fill the satellite's field of view ERR 
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cover the entire sky; therefore, it is more likely that the 
satellite sees the same sky conditions as the observer. 

Misclassification of nimbostratus in the 17 Jan and 28 Mar 
cases are due to warm infrared temperatures above the pre- 
established thresholds causing the cloud decision to be 
stratocumulus. Ground-observed stratocumulus is either 
classified as nimbostratus, if the infrared temperatures are too 
cold, clear/stratus, if the temperatures are too warm, or 
cumulus, if the texture test produces a high standard deviation 
(i.e. edge of a stratocumulus cloud located near a synoptic 
station). 

Two exceptions to the good performance by the stratiform 
class are altostratus and stratus. Both reports of altostratus 
are classified as stratus by the NPS model due to the satellite 
sensor's detection of warm infrared temperatures below the pre- 
established thresholds for altostratus (256 K, 17 Jan; 261 K, 20 
Mar). For both cases, the surface report is altocumulus 
translucidus at a single level (WMO middle cloud type 3, Table 
E. Due to the transparent nature of this type of middle cloud, 
surface radiation reaches the sensor causing the 
misclassification. 

There are also model-observation disagreements with stratus. 
In all five cases, this classification error can be attributed to 
infrared temperatures that are too cold; the temperatures exceed 
the thresholds and are specified as altostratus or nimbostratus. 
This error is an example of the limitations of objectively 
evaluating cloud detection. A difference in perspective between 


the ground observer (looks up) and the satellite (looks down) 
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causes this stratus misclassification. The ground observer is 
unable to detect the upper cloud layers that the satellite senses 
due to the opaque nature of the lowest cloud layer (Bunting and 
Hardy, 1984;  Tsonis, 1984); therefore, the surface report of 
stratus (lowest layer) and the NPS model cloud-typing decision of 
nimbostratus/altostratus do not agree. 

Misclassifications also occur with nonuniform or rough- 
textured clouds such as cumulus, cumulus congestus and 
cumulonimbus. The incorrect NPS model cloud type classifications 
are due to either low albedo values, warm infrared temperatures 
Gr Error in texture determination. 

Many of the surface reports of cumulus have visible albedo 
values less than the visible cloud/no cloud threshold (0.15). 
The small cumulus do not fill the satellite's field of view 
resulting in a classification of clear skies. The sensor is 
unable to resolve the individual cumulus clouds that the observer 
is able to see. Direct comparison between the surface 
observation and corresponding satellite pixel produced better 
results than the 7 x 7 matrix voting procedure (Table 4). When 
the field of view is smaller, the small cumulus elements are more 
likely to be resolved. This illustrates another difficulty in 
relating the ground-based observation to data collected from 
satellites. 

Many cumulus, cumulus congestus and both cumulonimbus reports 
also have infrared temperatures that are warmer than the pre- 
established thresholds due to the influence of surface radiation. 


Therefore, cumulus is classified as clear skies and cumulus 
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congestus and cumulonimbus are classified as either cumulus or 
clear skies. The majority of surface observations of cumulus 
congestus, analyzed as cumulus by the model, are ship reports. 
These ship reports are probably classified correctly as cumulus 
since cumulus congestus is a cloud type that is more likely to 
exist over land (more surface heating) than over water. 

Some surface reports of cumulus are classified as 
stratocumulus or stratus due to low standard deviation values 
signifying smooth textures. Changing the standard deviation 
criteria did not solve this problem since other correctly- 
classified stratocumulus and stratus clouds were then incorrectly 
categorized as cumulus. 

The disagreements between the NPS model output and the surface 
observations of cirrus are also due to extremely warm infrared 
temperatures. Many of the ground-observed cirrus-type clouds are 
EC TU S fibratus (WMO high cloud 1) or cirrostratus not 
progressively invading the sky (WMO high cloud 8), which are 
nearly transparent (Table 3). Optically thin cirrus are 
semitransparent to the upward directed radiation which allows the 
underlying area and clouds to contribute to the cloud-typing 
decision. All three correct cirrus classifications are of the 
type cirrostratus covering the entire sky (WMO high cloud 7) 
which fills the sensor's field of view; the cold infrared 
temperatures associated with this type of cirrus are well below 
the thresholds. 

An improvement in the cloud-typing statistics occurs when 
similar types of clouds are merged together such as the low 


clouds, cumulus, stratocumulus, stratus and cumulus congestus and 
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middle cloud type altostratus and nimbostratus. For example, in 
the 20 Mar case there are only two correct classifications in the 
direct surface-satellite comparison and 7 x 7 matrix (Table 4). 
Five out of twelve (42%) correct classifications are made when 
the low clouds are considered as one group (Table 5). The coarse 
resolution of the satellite data relative to the surface observer 
and the difference in viewing geometry limits the number of 
categories that can be distinguished by the NPS satellite 
algorithm (Bunting and Hardy, 1984). 

The majority of surface reports are multi-layered clouds. 
Since no separate distinction for multiple cloud layers is made 
by the NPS model, a special evaluation procedure is utilized. 
Most of the NPS model decisions are in agreement with the surface 
categorizations. For example, in the 13 May case, of the twenty- 
one surface reports of multi-layered clouds, there are nine 
belonging to the THIN category, four of type CL, three of type 
ML, one of type CM and CH and only three mismatches (Table 7). 

There are only two causes of mismatches between the satellite 
and conventional data: the edges of multiple cloud systems are 
located near synoptic stations which results in satellite 
decisions of clear skies or low clouds; and reports of low and 
high clouds are classified as middle cloud since the sensor 
detects a large dense cloud mass and is unable to distinguish 
between the different layers. The NPS model does a fair job 
portraying the multi-layered cloud systems. 

The objective precipitation statistics are presented in Table 


e The NPS model is able to delineate most precipitation areas 
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since the majority of surface observations of precipitation are 
identified by the model. The total percentage correct value for 
all reports of precipitation is 49% with a range of values from 
41% in the 17 Jan Case to 100% in the 13 May case. 

The discrepancies between the model and ground truth data are 
caused by infrared temperatures that are warmer than the 
precipitation/no precipitation threshold. For example, in the 17 
Jan case, several observers report precipitation with 
stratocumulus clouds. The NPS model classifies the stratocumulus 
correctly, but does not detect the precipitation since the 
infrared temperatures are warmer than the precipitation/no 
precipitation threshold of Zag 5 Ke. 

Ninety-four percent of all surface reports Of no 
precipitation are in agreement with the satellite analysis (Table 
7). The most successful result is 97% in the 13 May case and the 
least successful percentage correct value (90%) occurs in the 17 
Jan study. The satellite-surface disagreement can be attributed 
to Visible and infrared values exceeding their respective 
thresholds. 

In general, the objective evaluation results of cloud type and 
precipitation identification are encouraging. The following 
chapter presents the results of the four independent 


investigations used to test the new thresholds. 
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VEITI: STATISTICAL RESULTS OF THE INDEPENDENT TESTEO ORR 


The statistical results of four independent cases (20 Feb , 
27 Mar , 16 Apr and 13 May) used to test the hypothesized 
thresholds are presented in Tables 8-11. Each new study utilizes 
the corresponding tuned thresholds of the previous 
investigations: The 20 Feb case uses 20 Mar thresholds; the 27 
Mar case uses the Mar 28 values; the 16 Apr case uses 11 Apr 
thresholds; and the 15 May case uses the 13 May values. 

The results of the objective evaluation of the independent 
studies are similar to those of the initial cases. The most 
successful cloud type classification (78% correct for all cases) 
clear skies category. Percentage correct values range from 63% 
(20 Feb case) to 100% (16 Apr and 15 May case) with better or at 
least as good statistics obtained with the 7 x 7 voting scheme 
(Table 8). The low percentage of the 20 Feb study is related 
to clear skies over snow-covered ground that the NPS model 
identifies as low clouds. 

The best single cloud type classifications are for stratiform 
types nimbostratus and stratocumulus. Sixty percent of all 
nimbostratus reports are classified correctly by the model; 
percentage correct values range from 25% (15 May case) to 100% 
(20 Feb case). For stratocumulus the range of values is 29% 
(27 Mar case) to 50% (16 Apr case) with a total value of 32% for 
all cases. The statistical outcomes from the 7 x 7 voting 
procedure are better or as good as the direct pixel- to- station 


comparison for these uniform-textured clouds. Most of the 
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disagreement between the model and surface data in the 
stratocumulus classification is due to the difference in 
perspective between the ground and satellite. The surface 
observer views the low cloud stratocumulus while the satellite 
senses multiple layers. The ground observer can not identify the 
upper cloud layers. 

As in the initial cases, the model has difficulty specifying 
the other stratiform cloud types stratus and altostratus. Only 
one of the nine reports of stratus is classified correctly by the 
model. Most of the stratus, analyzed incorrectly as clear skies, 
are located on the edges of synoptic stations. The one report of 
altostratus is analyzed as nimbostratus by the model; the 
infrared temperatures are colder than the threshold. 

Nonuniform-textured cumulus, cirrus, cumulus congestus and 
cumulonimbus clouds are not easily identified by the model in the 
test cases. Three out of thirty-one reports of cirrus and two 
out of eight reports of cumulus are classified correctly. 
Cumulus congestus and cumulonimbus are never accurately depicted. 

The same classification errors for these  rough-textured 
clouds are observed in the test cases as in the initial studies. 
The cirrus is semi-transparent and the Cumulus do not fill the 
sensor's field of view; surface characteristics influence the 
sensor's interpretation of the cloud's temperature and albedo 
causing the analysis error. The five surface reports of 
cumulonimbus and thirteen reports of cumulus congestus are 
misclassified as cumulus congestus and cumulus, respectively, 
when the infrared temperatures are warmer than the predetermined 


thresholds. Cumulonimbus is analyzed as nimbostratus when the 
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visible albedos are below the 0.88 visual threshold separating 
these two cloud types. 

The majority of surface reports are of multiple cloud layers 
and are in agreement with the subjective evaluation of the 
satellite data. Out of eighty-six multi-layered cloud reports, 
only five are labelled mismatch--the satellite and surface data 
are in disagreement (Table 10). The majority of surface reports 
of multiple cloud layers correspond to a satellite matrix that 
contains classifications of cloud types at different levels 
(MULTE). 


The precipitation statistics for the independent cases are 


presented in Table 11. Elghty-three percent of the surface 
reports of precipitation are recognized by the model. The no- 
precipitation results are also reasonable. Ninety-three percent 


of all no-precipitation cases are correct with values ranging 
from 84% (27 Mar case) to 98% (20 Feb case). As in the initial 
cases, the only surface-satellite disagreement of precipitation 
occurrence in the test cases is due to visible and infrared 
values above or below the precipitation/no precipitation 
thresholds. 

The cloud and precipitation statistical results of the 
independent test cases resemble the outcomes of the initial 
investigations. Similar classification successes and failures 
are noted. The tuned thresholds from the initial studies were 
successful fairly well for the test cases. Generalizing the NPS 
model's algorithm by altering the infrared thresholds to account 


for the observed monthly variation in surface radiation produces 


46 


reasonable cloud and precipitation analyses. The successes and 
limitations of the model as well as recommendations for future 


research are described in the following chapter. 
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IX. CONCLUSIONS AND RECOMMENDATIONS 


Directly integrating satellite and surface data is a 
difficult task. Satellite imagery provides important cloud and 
precipitation information from above while conventional surface 
data describes weather features from below. This thesis 
illustrated that using satellite and surface data as 
complimentary forecasting tools produces better analyses of 
mesoscale and subsynoptic weather features. The performance of 
an automated cloud and precipitation model (Wash et al., 1985) 
that uses GOES visible and infrared data to provide high 
resolution analyses was evaluated for several 1985 winter and 
spring days. The infrared thresholds used in the algorithm were 
modified for each case due to the observed temperature 
variability from January to May. 

The NPS model succeeds in analyzing most overcast and 
clear skies but underestimates broken and scattered situations. A 
visible albedo of 0.15 seems to be a representative cloud/ no 
cloud threshold. Most surface reports of clear skies are in 
agreement with the model except clear skies over snow-covered 
ground which is analyzed as low cloud since the visible and 
infrared thresholds are exceeded. 

The cloud type analyses are reasonable. Accuracy increases 
when all low level clouds or middle level clouds are considered 
as one cloud type. The majority of large uniform-textured 
clouds such as nimbostratus and stratocumulus are classified 


correctly by the model. The 7 x 7 voting scheme is superior to 
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the direct pixel-station comparison method when evaluating these 
stratiform clouds. The majority of surface reports of multi- 
layered clouds are also described correctly by the model. 

Some stratiform and cumuliform clouds are misclassified due 
to an error in texture decision. Using more than four albedos 
when computing the standard deviations used in the texture test 
may reduce this error. 

Although some large-sized cumuliform clouds are identified 
by the model, small cumulus are not distinguished by the model; 
they are classified as clear skies. Since they fail to fill 
the instrument's field of view, surface radiation reaches the 
sensor causing a misclassification. Direct pixel-station 
comparison statistics are better than the 7 x 7 voting 
procedure results with cumuliform clouds. 

Some dense cirrus are analyzed correctly by the model but 
semitransparent cirrus are not detected: by the model and are 
identified as clear skies or low clouds. Again, surface 
radiation influences the cloud-typing decision. 

Most cumulus congestus clouds are analyzed as cumulus since 
the infrared temperatures are warmer than the threshold used to 
separate cumulus congestus from cumulus. It is probable that the 
majority of these reports of cumulus congestus are analyzed 
correctly as cumulus since they are ship observations; cumulus 
congestus is more likely to exist over land than water. 

Cumulonimbus are classified incorrectly as cumulus, cumulus 
congestus or nimbostratus since the infrared temperatures and 
Visible albedos are too warm and/or too dim, respectively, to 


result in the cumulonimbus classification. 
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The cloud-top temperature/height analyses are generally 
accurate. The model is able to depict the wide range of 
temperatures/heights associated with the various cloud types of 
a frontal Zone. 

The NPS model produces fair analyses of precipitation 
distribution. Some areas of precipitation are not identified by 
the model since the visible and/or infrared values are below 
the precipitation/no precipitation thresholds. Some 
nonprecipitating regions are identified as precipitación areas 
since the visible and/or infrared values are above the 
precipitation/no precipitation thresholds. 

Some error in estimating precipitation intensity exists. 
Many times the intensity is underestimated. This is a 
difficult quantity to measure using digital satellite data. 

Discrepancies between the surface and satellite data may 
be caused by the difference in perspective between the ground and 
satellite, alignment of the cloud to the surface stations and 
misleading WMO code descriptions. Many surface reports of 
dense low clouds are classified as nimbostratus (multi-layered) 
by the model. The surface observer cannot see beyond the 
lowest cloud layer. The satellite sees the top of the cloud 
structure. An observer views clouds subjectively at different 
angles for variable areal regions. A satellite interprets 
objectively the infrared and visible information at a single 
point: Therefore, problems occur with direct comparison of 


satellite pixel information and surface observations. 


50 


Many times clouds are located on edges of synoptic stations 
which leads to surface-satellite disagreement. Navigation error 
may also contribute to this misalignment. 

Some of the WMO code descriptions combine several cloud 
types under one category (i.e. middle cloud WMO code 7) which 
contributes to subjective cloud interpretation differences. 

Recommendations for future research include the following: 
(1) differentiating between snow cover and clouds by 
implementing a spectral and spatial scheme (Tsonis, 1984) that 
separates classes of cloud and surface characteristics; (2) 
using more than four points when computing cloud amount and 
standard deviations used in the texture test which may reduce 
the broken and scattered and stratiform/cumuliform classification 
errors; (3) optimizing the model so that coarser visible and 
infrared resolution can be used to produce cloud and 
precipitation analyses; (4) further testing of the threshold 
values using GOES data from different geographical regions; 
NS utilizing digitized GOES-TAP data as input compared 


to stretched digital VISSR data with the NPS model. 


Sai 


Code Figure 


O 


APPENDIX A 


TABLES 


TABLE 1 


WMO Low Cloud Types 


Specification 


No Stratocumulus, Stratus, Cumubuss 
or Cumulonimbus 


Cumulus humulis or Cumulus fractus other 
than that of bad weather, or both 


Cumulus mediocris or congestus, with or 
without Cumulus of specis fractus or 
humulis or Stratocumulus, all having their 
bases at the same level 


Cumulonimbus calvus, with or without 
Cumulus, Stratocumulus or Stratus 


Stratocumulus cumulogenitus 


Stratocumulus other than Stratocumulus 
cumulogenitus 


Stratus nebulosus or Stratus fractus other 
than that of bad weather, or both 


Stratus fractus or Cumulus fractus of 
bad weather, or both (pannus), usually 
below Altostratus or Nimbostratus 


Cumulus and Stratocumulus other than 
Stratocumulus cumulogenitus, with bases at 
different levels 


Cumulonimbus capillatus (often with an anvil) 


with or without Cumulonimbus calvus, Cumulus, 
Stratocumulus, Stratus or pannus 
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TABLE 2 


WMO Middle Cloud Types 


Specification 
No Altocumulus, Altostratus or Nimbostratus 
Altostratus translucidus 
Altostratus opacus or Nimbostratus 
Altocumulus translucidus at a single level 


Patches (often lenticular) of Altocumulus 
translucidus, continually changing and 
occurring at one or more levels 


Altocumulus translucidus in bands, or one 

or more layers of Altocumulus translucidus 
or opacus, progressively invading the sky; 
these Altocumulus clouds generally thicken 
as a whole 


Altocumulus cumulogenitus (or cumulonimbo- 
genitus) 


Altocumulus translucidus or opacus in two or 
more layers, or Altocumulus opacus in a 
single layer, not progressively invading 

the sky, or Altocumulus with 

Altostratus or Nimbostratus 


Altocumulus castellanus or floccus 


Altocumulus of a chaotic sky, generally at 
several levels 
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Code Figure 


0 


i 


TABLE 3 


WMO High Cloud Types 


Specification 
No Cirrus, Cirrocumulus or Cirrostratus 


Cirrus fibratus, sometimes uncinus, not 
progressively invading the sky 


Cirrus spissatus, in patches or entangled 
sheaves, which usually do not increase and 
sometimes seem to be the remains of the 
upper part of a Cumulonimbus; or 
castellanus or floccus 


Cirrus spissatus cumulonimbogenitus 


Cirrus uncinus or fibratus, or both, 
progressively invading the sky; they generally 
thicken as a whole 


Cirrus (often in bands)and Cirrostratus, or 
Cirrostratus alone, progressively invading 
the sky; they generally thicken as a whole, 
but the continuous veil extends more than 
45 degrees above the horizon 


Cirrus (often in bands)and Cirrostratus, or 
Cirrostratus alone, progressively invading 
the sky; they generally thicken as a whole, 
but the continuous veil extends more than 
45 degrees above the horizon, without the 
sky being totally covered 


Cirrostratus covering the whole sky 


Cirrostratus not progressively invading the 
sky and not entirely covering it 


Cirrocumulus alone, or Cirrocumulus predomi- 
nant among the high clouds 
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Direct Comparison 


err 


Ci 
As 
SC 
SC 
Ns 
Cu 
ec 
ED 


cir 
aan 
AS 
st 
SE 
Ns 
Cu 
CE 
Cb 


TABLE 4 


Specific Cloud Type Statistics for Initial Cases 


17 Jan 


4/4 (100%) 
no report 
0/1 
2/8 (25%) 
10/25 (403) 
6/21 (29%) 
0/4 
0/1 
no report 


28 


9/10 (90%) 
1/1 (1003) 
no report 

0/2 

8/20 (40%) 
4/7 (57%) 

3/5 (60%) 

0/2 

no report 


ES 


Clr 24/25 (96%) 


Ci 
AS 
SE 
sc 
Ns 
Cu 
CC 
ep 


Wy Gall 7 3) 
no report 
071 
0/1 
no report 
1/37333) 
079 
0/1 


4/4 (100%) 
no report 
0/1 

1/8 (133) 

14/25 (56%) 
8/21 (38%) 
0/4 

0/1 

no report 


Mar 


9/10 (90%) 
1/1 (1003) 
no report 
0/2 

117/20 (553) 
4/7 (57%) 

2/5 (40%) 

07/2 

no report 


May 


7 x 7 Vote 


25/25 (100%) 


AGE) 
no report 
0/1 
0/1 
no report 
1/3 (33%) 
0/9 
0/1 
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Direct Comparison 


20 Mar 


19/21 (90%) 
2/11 (183) 
0/1 
no report 
1/5 (203) 
no report 
1/4 (25%) 
0/3 
no report 


11 Apr 


als) 

0/5 

no report 
0/1 

5/14 (363) 
no report 
0/4 

0/1 

no report 


7 x 7 Vote 


19/21 (90%) 
P I1 

0/1 

no report 

1/5 (20%) 

no report 
1/4 (253) 

0/3 

no report 


6/7 (86%) 

0/5 

no report 

0/1 

6/14 (433) 
no report 

0/4 

0/1 

no report 


Total for all cases 


61/67 ( 
A23 
0/2 
DARI 

24/65 

10/28 
5/20 
0/16 
0/1 


63/67 
3223 
0/2 
11,20, 

32/65 (4 

I/28 (4 
4/20 (2 
0/16 
0/1 


TABLE 5 


Cloud Group Statistics for Initial Cases 


Direct Comparison 


Cir 
CH 
CM 
CLG 
CB 


Cir 
CH 
CM 


CB 


Cir 


CM 
CE 
CB 


4/4 (100%) 
no report 
11/22 (50%) 
24/38 (63%) 
no report 


19/21 (903) 
17/11 
0/1 
5/12 (423) 
no report 


9/10 (90%) 

1/1 (1005) 

B DAE 
19/29 (663) 
no report 


S070 67 1$) 
0/5 
no report 
10/20 (50%) 
no report 


24/25 (96%) 
1/6 17s) 

no report 
4/14 (29%) 
07 1 


61/67 (91%) 
2223 B RS) 
16/30 (532) 
62/113 (552) 

0/1 


17 


20 


28 


rI 


ES 


Total for all cases 


Jan 


Mar 


Mar 


Apr 


May 
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7 x 7 Vote 


4/4 (100%) 
no report 
10/22 (45%) 
25/38 (66%) 
no report 


19/21. (90%) 
le: 
O7 a 
4/12 (333) 
no report 


9/10 (90%) 

1/1 (1005) 

4/7 (57%) 
19/29 (66%) 
no report 


5/7 ia) 
O75 
no report 
10/20 (50%) 
no report 


25/25 (1003) 
AS RS 

no report 
5/14 (363) 
0/1 


63/67 (94%) 
3/23 (Nati 
14/30 (47%) 
63/113 (56%) 

0/1 


TABLE 6 


Multi-layered Cloud Statistics for Initial Cases 


Percentage satellite- 
surface agreement 


Percentage satellite- 
surface mismatched 
Percentage satellite- 
surface agreement 
Percentage satellite- 
surface mismatched 
Percentage satellite- 


surface agreement 


Percentage satellite- 
surface mismatched 


Percentage satellite- 
surface agreement 


Percentage satellite- 
surface mismatched 
Percentage satellite- 
surface agreement 
Percentage satellite- 
surface mismatched 
Percentage satellite- 


surface agreement 


Percentage satellite- 
surface mismatched 


15/17 


27/087 


18/20 


2/20 


17/19 


2/19 


30/33 


37/33 


18/21 


Bye 


17 Jan 
(88%) 


(123) 


20 Mar 
(90%) 


(10%) 


28 Mar 
(89%) 


(113) 


Apr 11 
(91%) 


(9%) 


13 May 
(863) 


(143) 


Total” tor 
98/110 (89%) 


12/110 (T13) 


27 


THIN 
CH 

CM 
CL 
MULTI 


THIN 
CH 

CM 
CE 
MULTI 


THIN 
CH 

CM 
CL 
MULTI 


THIN 
CH 

CM 
CL 
MULTI 


all cases 


THIN 
CH 
CM 
en 
MULTI 


DAT 
By 
M17 
2417 
QUE 


3/20 
3/20 
0/20 
3/20 
9/20 


6/19 
2/19 
1/19 
4/19 
4/19 


6/33 
8/33 
1/38 
4/33 
11/33 


9/21 
4/21 
0201 
10020 
3/21 


26/110 
20/110 

4/110 
14/110 
34/110 


> 


TABLES? 


Precipitation Statistics for Initial Cases 


17 Jan 
PRECIPITATION 9/22 (413) 


NO PRECIPITATION 54/60 (903) 


20 Mar 
PRECIPITATION no rain 


NO PRECIPITATION 61/65 (943) 


28M 
PRECIPITATION 7/11 (643) 


NO PRECIPITATION 50755 009 9 


ll Apr 
PRECIPITATION 2/31 (61%) 


NO PRECIPITATION 64/68 (943) 


13 May 
PRECIPITATION 1/1 (100%) 


NO PRECIPITATION 65/67 (97%) 


Total for all cases 


PRECIPITATION 18/37 (49%) 


NO PRECIPITATION 295/315 (943) 
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Direct Comparison 


elr 
Ci 
AS 
su 
Sc 
Ns 
Cu 
ec 
CD 


cir 
Ci 
AS 
st 
Sc 
Ns 
Cu 
ec 
ED 


Cile 
Ci 
AS 
St 
Sc 
Ns 
Cu 
eC 
Cb 


eur 
ET 
AS 
St 
SC 
Ns 
Cu 
GE 
CD 


err 
Ci 
As 
St 
SE 
Ns 
Cu 
CC 
Cb 


TABLE 8 


20 Feb 
41/65 (63%) 
0/5 
no report 
1/6 (17%) 
3/10 (30%) 
2/2 (100%) 
no report 
no report 
no report 
27 Mar 
16/17 (94%) 
3/20 ml 5) 
no report 
0/2 
4/14 (293) 
5/6 (835) 
0/3 
0/1 
no report 
16 Apr 
19/19 (100%) 
no report 
no report 
0/1 
2/4 (50%) 
1/3 (33%) 
2/5 (40%) 
0/6 
0/4 
15 May 
12/12 (100%) 
0/6 
0/1 
no report 
3/9 (333) 
1/4 (25%) 
no report 
0/6 
o 
Total for all cases 
88/113 (78%) 
3/31 (10%) 
0/1 
1/9 (113) 
12/37 (323) 
9/15 (60%) 
2/0 (25%) 
0/13 
0/5 
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Specific Cloud Type Statistics for Independent Test Cases 


ux 7 Vote 


45/65 (69%) 
075 
no report 
0/6 
3/10 (30%) 
2/2 (1003) 
no report 
no report 
no report 


17/17 (100%) 
3/20 (15%) 

no report 
0/2 
4/14 (293) 
5/6 (83%) 
1/3 (335) 
O/I 

no report 


19/19 (100%) 

no report 

no report 
0/1 
2/4 
1/3 
179 
0/6 
0/4 


50 
33 
20 


A PA a S 


o? o9 oe 


12/12 (1003) 
0/6 
0/1 

no report 
5/9 (563) 
1/4 (25%) 

no report 
0/6 

01 


93/113 (827) 
2/31 (103) 
0/1 
0/9 

14/27 391) 
8/15 (53%) 
2/8 (25%) 
0/13 
0/5 


TABLE 9 


Cloud Group Statistics for Independent Test Cases 


Direct Comparison 


Cir 
CH 
CM 


CB 


CLE 
CH 
CM 


CB 


7 Xx lute 


20 Feb 
41/65 (63%) 45/65 (693) 
0/5 0/5 
2/2 (1005) 1/2 (50%) 
6/16 (383) 5/16 31%) 
no report no report 
27 Mar 
17/17 (10037) 17/17 (1005) 
3/20 (155) 3/20 (15%) 
5/6 (833) 5/6 (833%) 
5/2070 (2537) 5/20 (2957) 
NO report no report 
16 Apr 
19/19 (1003) 19/19 (1003) 
no report no report 
2/3. (67%) 2/3 (67%) 
l0 516 70635) 10/16 (63%) 
0/4 0/4 
15 May 


12712 (1002) 
0/6 

4/5 (80%) 
10/15 (67%) 
0/1 


89/113 (79%) 
3/30 (102) 

13/16 (81%) 

31/67 (46%) 

0/5 


12/12 (1002) 
0/6 

4/5 (80%) 
10/15 (67%) 
0/1 


Total for all cases 
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93/11 3825 
3/31 (10%) 

12/16 (75%) 

30/67 (45%) 

0/5 


TABLE 10 


Multi-layered Cloud Statistics for Independent Test Cases 


Percentage satellite- 
surface agreement 


Percentage satellite- 
surface mismatched 


Percentage satellite- 
surface agreement 


Percentage satellite- 
surface mismatched 


Percentage satellite- 
surface agreement 


Percentage satellite- 
surface mismatched 


Percentage satellite- 
surface agreement 


Percentage satellite- 
surface mismatched 


Percentage satellite- 
surface agreement 


Percentage satellite- 
surface mismatched 


16/17 


P207 


25/28 


3/28 


290/21 


eu 


20/20 


0/20 


20 Feb 


(94%) 


27 Mar 


(89%) 


(113) 
16 Apr 


(95%) 


15 May 


(100%) 


THIN 


CM 
CE 
MULTI 


THIN 
CH 
CM 


MULTI 


THIN 
CH 

CM 
CL 
MULTI 


Total for all cases 


81/86 


5/86 


(94%) 


(6%) 


6l 


THIN 
CH 

CM 
CL 
MULTI 


5/17 
0/17 
i717 
252p 
3/17 


6/28 
37008 
0/28 
0/28 
16/28 


272] 
0/21 
17/21 
4/21 
13/21 


4/20 
102219 
1/20 
9/20 
5/20 


17/86 
4/86 
3/86 

20/86 

37/86 


TABLE 11 


Precipitation Statistics for Independent Test Cases 


20 Feb 
PRECIPITATION 2/3 (673) 
NO PRECIPITATION 100/102 (98%) 

27 Mar 
PRECIPITATION 11/11 "CLOS 
NO PRECIPITATION 67/80 (84%) 

16 Apr 
PRECIPITATION 4/6 (67%) 
NO PRECIPITATION 55/57 (963) 

15 May 
PRECIPITATION 3/4 (753) 
NO PRECIPITATION 52/55 (961) 


Total for all cases 


PRECIPITATION 20/24 (833) 
NO PRECIPITATION 274/294 (933) 


62 


APPENDIX B 


FIGURES 
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Figure 1. Two-dimensional Cloud-typing Graph Using GOES 


Infrared and Visible Satellite Digital Data for 
i7 Jan 85 
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Two-dimensional Cloud-typing Graph Using GOES 
Infrared and Visible Satellite Digital Data for 


20 Mar 85 


Figure 2. 


64 


165 


Ns/MI Co 
Ci 
254 


AS 
e K Cu Cong 


(temperature) 


K 
286 


St/Fog 


291 





AS 1.0 


VIS 


(albedo) 


Figure 3. Two-dimensional Cloud-typing Graph Using GOES 
Infrared and Visible Satellite Digital Data 
for 28 Mar 85 
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Figure 4. Two-dimensional Cloud-typing Graph Using GOES 
Infrared and Visible Satellite Digital Data 


for 11 Apr 85 
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Figure 5. Two-dimensional Cloud-typing Graph Using GOES 
Infrared and Visible Satellite Digital Data 


for 13 May 85 
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Figure 6. Two-dimensional Cloud-typing Graph Using GOES 
Infrared and Visible Satellite Digital Data 
for Summer 83 
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Observed Stratocumulus 


Se tor Se tor Se tor 
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Figure 7. Scattergram of Visible and Infrared Raw Data Counts 
Corresponding to Surface Reports of Stratocumulus 
(Lower Infrared Count for Warmer Temperature) 
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Observed Cumulus and Stratus 


Cu and St Cu and St Cu and St 
for 20 MAR for 28 MAR for 11 APR 
T X 





Figure 8. Scattergram of Visible and Infrared Raw Data Counts 
Corresponding to Surface Reports of Cumulus and 
Stratus (Lower Infrared Count for Warmer Temperature) 
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Observed Multi-layered Clouds 


ML tor ML. for ML. for ML for ML for 
17 JAN 20 MAR 28 MAR t1 APR 13 MAY 
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Figure 9. Scattergram of Visible and Infrared Raw Data Counts 
Corresponding to Surface Reports of Multi-layered 
Clouds (Lower Infrared Count for Warmer Temperature) 
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Figure 10. Two-dimensional Frequency Graph Using GOES Infrared 
and Visible Satellite Digital Data for 17 Jan 85. 
The Number in the Darkened Regions Represents the 
Average Frequency of Occurrence (Higher Frequency 
in Darker Areas) of Each Infrared and Visible (Scaled 
to Infrared Data by Mulitplying by Four) Count 
Combination That Occurs Within the Image 





Figure ll. Same as Fig. 10 except for 20 Mar 85 
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Figure 12. Same as Fig. 10 except for 13 May 85 
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Figure 14. GOES Infrared Imagery for 1800 GMT 17 Jan 85 
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Cloud Amount Analysis for 1800 GMT 17 Jan 85. Four 
Colors Are Used To Distinguish Between Cloud Amounts: 
NOS Cloud Cover (OVE Blue; 75% Cloud Cover 
CER O Ci bark Green; 50% Cloud Cover (BRO), 

Light Green; 25% Cloud Cover (SCA), Red 
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eeud i peana lysis for 1800 GMT 17 can 85. Eight 
Colors are Used TO Tllustrate Eight Different Cloud Type 
Cumulonimbus ch) Dark Blue: Mranbostratus (Ns), Blue; 
nulle congestus Ace) Light Blue; Cumulus (Cu), Dark 


Green; Stratocumulus (Sc), Light Green; Stratus (St), 
Mello Na leostratus (As), Orange; Cirrus (Ci), Red 


Figure l7; 





Cloud-top Temperature Analysis for 1800 GMT 17 Jan. 

Nine Colors Are Used To Illustrate 10 K Intervals of 
Cloud-top Temperatures: 210-219 K, Dark Blue; 220-2900 
Blue; 230-239 K, Light Blue; 240-249 K, Dark Green; 
250-259 K Green; 260-269 K, Light Green; 270-279 K, Valim 
280-289 K, Red; 290-300 Kadsursacei ca 
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Figure 18. Automated Radar Summary Chart for 17 JAN 85 
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Figure 19. Cloud-top Height Analysis for 1800 17 Jan 85. Nine 
Colors Are Used To Illustrate 100 mb Cloud-top Heres 
Intervals: 100-199 mb, Dark Blue; 200-299 mb, Blue; 
300-399 mb, Light Blue; 400-499 mb, Dark Green: 
500-599 mb, Green; 600-699 mb, Light Green; 700-799 mb, 
Yellow; 800-899 mb, Red; 900-1000 mb; Dark Red 


ze 

Figure 20- COS S Teen ede for 1800 GNP 
17 Jan 85. Three Colors Are Used To Distinguish 
Between Precipitation Intensities: Light, Red; 
Moderate, Green; Heavy, Blue 
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Synoptic Surface Station Observations of WMO Cloud 
Types and Present Weather Over GOES Visible Imagery 

for TSC0 CHIP IR Jjançãs: The Number In the Upper Lert 
Corner Is the WMO Low Cloud Type; the Number In the 
Lower Left Corner Is the WMO Middle Cloud Type; and the 
Number In the Upper Right Corner Is the WMO High 

Cloud Type. Standard WMO Weather Abbreviations 

Are Found In the Lower Right Corner. 
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Figure 23. GOES Visible Imagery for 20 Mar 85 
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Figure 24. GOES Infrared Imagery for 20 Mar 85 
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Figure 25. Same as Fig. 15 except for 11 Apr 85 
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Dro; GOES eintraned Imagery for 11 Abr 85 
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Figure 28. Snow Cover In Inches for 1200 GMT 10 Apr 85 
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Figure 29. Same as Fig. 16 except for 28 Mar 85 
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Figure 30. GOES Visible rmadery dor 28 Tm RE 
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Figure 31. GOES Infrared Imagery for 28 Markso 
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Figure 32. Same as Fig. 16 except for ll Apr 85 


87 


Pa s LT uc 


4 3 


; 
Af 


A 





Figure 


sd 
Ei 
» 

E 
ES 
EF 
re: 
E 
t 
Y 
, 





Figure 34. GOES Infrared Imagery for 13 May 85 
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Figure 35. Same as Fig. 16 except for 13 May 85 
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Figure 36. Same as Fig. 17 except for 28 Mar 85 
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Figure 37» Same as Fig. 20 exccoron oO Marias 
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Figure 38. Automated Radar Summary Chart for 20 Mar 85 
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Figure 39. Two-dimensional Precipitation Intensity Estimate Graph 
Using GOES Infrared and Visible Satellite Digital Data 
fOr Zonar do 
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Figure 40. Same as Fig. 20 except for ll Apr 85 
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Figure 41. Automated Radar Summary Chart for 11 Apr 85 
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